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Abstract 

It was recently shown that entanglement in quantum systems being in a non-equilibrium state can 
appear at much higher temperatures than in an equilibrium state. However, any system is subject 
to the natural relaxation process establishing equilibrium. The work deals with the numerical study 
of entanglement dynamics in a dipolar coupled spin-1/2 system under the transition from a non- 
equilibrium state to an equilibrium state. The spin system is characterized by a two-temperature 
density matrix, and the process of the establishment of equilibrium is in the equalization of these 
temperatures. The method of the non-equilibrium statistical operator is used to describe the 
evolution of the system. The process of establishing an equilibrium state in the homonuclear spin 
systems at low temperature was first considered. It was shown that the time dependences of 
the inverse temperatures of the spin subsystems are given by a solution of non-linear equations 
in contrast to the linear equations in the well-known high temperature approximation. It was 
first studied the entanglement dynamics during the equilibrium state establishment and found 
that, during establishing equilibrium, the concurrence changes non-monotonically with time and 
temperatures. Entanglement fades long before equilibrium is established in the system. It was 
shown that the entanglement dynamics depends strongly on the ratio of the Zeeman energy to the 
dipolar energy. At a high ratio, the concurrence in the system decreases quickly for time about 
100 fis while, at a low ratio, establishment of equilibrium and fading entanglement take prolonged 
time up to 1 ms. 
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I. INTRODUCTION 



Entanglement is a term used in quantum theory to describe the fact that the particles 



can be correlated with^each other regardless of how far they are from each other 
Many current studies 



focus on how to harness the potential of entanglement in the 
development of quantum cryptography, quantum communication ^ , quantum computation 



and quantum metrology |6 



The growth in interest in entanglement applications has 



stimulated intensive qualitative and quantitative research of entang! 



systems and in both equilibrium 



pseudo-pure state |18|-|20| initial state 



3, QBE 



ement in various physical 



12 [and non-equilibrium 



13 



-17 



22|, including the 



It was shown that quantum systems in thermal equilibrium became entangled, if the in- 
teraction energy between them is larger than thermal energy due to their coupling to the 
environment. This condition leads to that entanglement appears only at very low tempera- 
ture. 

For example, entanglement in dipolar coupled spin system in equilibrium state is achieved 
by the application of a low external magnetic field in which the Zeeman interaction energy is 
the order of or even less than the dipolar interaction energy, ojq < uod 
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ll| . For a proton 



spin systems, the dipolar energy is Ud = 15 kHz. That results in the temperature of the 

-|l2|. Here is the Boltzmann constant. 



entanglement appearance of T ~ = 0.7 fiK 

Recently, it was proved that this relationship between temperature and entanglement 
is not justified for systems that are not in thermal equilibrium 23|]. In contrast to the 
equilibrium it was shown that entanglement in quantum systems being in a non- 



23 



24( 1 . Investigations of spin 



equilibrium state can appear at much higher temperatures 
systems demonstrate this statement: e.g. a spin system, which is in a separated state at 
equilibrium, is irradiated by a radiofrequency field and goes out from the equilibrium, and 



ISHlTl, 



22|. 



entanglement appears in the system simultaneously 

The present paper deals with the opposite case when the system being initially in a non- 
equilibrium state evolves to an equilibrium state. This entanglement dynamics during the 
natural process of establishing equilibrium is an actual problem which, according to our 
knowledge, does not studied yet. We study entanglement dynamics in a dipolar coupled 
spin-1/2 system under the transition from non-equilibrium state to equilibrium state. We 
consider a spin system, which can be divided to two thermodynamic subsystems which are 
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in the equilibrium state and characterized by two different temperatures (so-called a quasi- 



equihbrium state 



25 



26|). The process of establishing equilibrium in the whole system can 



be represented as equalization of their temperatures. Such quasi-equilibrium states were 



prepared experimentally using various methods |25N27l|. These states can be considered as 
the last stage in relaxation process from non-equilibrium state to equilibrium one. Existence 
of such states can be explained by using the fundamental principles of non-equilibrium 



statistical mechanics first formulated by Bogolyubov [28[. In non-equilibrium statistical 
mechanics, a fundamental role is played hierarchy of the relaxation times and reducing in 
the number of parameters requested for the description of a non-equilibrium process. As 
shown by Bogolyubov, after the lapse of a small interval of time, of the order of the duration 
of a collision, all higher distribution functions are completely determined by a single-particle 
function, and kinetic equations can be found for this stage . This means it is possible to 
use an abbreviated description of a non-equilibrium state in which the complete distribution 
function depends on few parameters. 

Bogolyubov's idea about hierarchy of the relaxation times has been widely used for de- 



scription of relaxation processes in spin systems 251. 129l-l31|. In particular, if the correlation 



times are much shorter than the time T2, the lifetime of the free precession signal, and the 
time T2 is shorter then the characteristic time of an investigated process, then the equilib- 
rium state is established in parts of the spin system before the whole system reaches the 
total equilibrium. 

The structure of the paper is as follows: in the next section, we describe the Hamiltonian 



for a spin system in an external field at the equi 



25 



ibrium and quasi-equilibrium states and 



31| . According to this theory in a dipolar 



principale of the Provotorov saturation theory 
coupled spin system in a quasi-equilibrium state, two spin temperatures and Td can be 
assigned to the Zeeman and dipolar thermodynamic subsystems 
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31| . In section 3 estab- 



lishment of the equilibrium state is analyzed. Then we consider dynamics of the pairwise 
entanglement in a quasi-equilibrium state which is characterized by these two temperatures. 
Discussion of the results is given in the final section. Numerical calculations of evolution of 
the spin temperature and concurrence Cmn between the m-th and n-th spins are performed 
using the software based on the MatLab package. To demonstrate dependence of the en- 
tanglement dynamics on the relative spin position we will consider three four-spin systems: 
linear chain, ring, and rectangle. The choice of these systems is due to the following reasons: 
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i) these structures exist in real molecules and used in experiments and ii) as dipole-dipole 
interactions decrease inversely to the third power of the distance between the spins, the 
effect of the system geometry is manifested more clearly at a small number of spins. 



II. SPIN SYSTEM IN QUASI-EQUILIBRIUM STATE 

We consider the system of identical spin- 1/2 particles with the dipole-dipole interaction 
(DDI) in an external magnetic field Hq along the z - axis. In the equilibrium state the 
system is characterized by a density matrix 

exp{-/3H) 
P= Z{P) ' 

with a single parameter P = denotes the inverse spin temperature, T is the spin tem- 
perature, and Z = Tr {exp (—/?//)} is the partition function. The Hamiltonian consists 
of the two parts, H = Hz + Hdd, the Zeeman part 

H, = -^H,Y,I^ (2) 
j 

and the dipolar part 

Hdd = -l^fi >J ( -B- ( ^J^Jk ) ( IkTjk ) - 1 • (3) 




Here 7 is the gyromagnetic ratio, Ij is the angular momentum operator of the j-th spin (j = 
1,2...N), fjk is the radius vector from the j-th to k-th spins with the spherical coordinates 
Tj-fc, 9jk and ipjk ■ Oj^ is the angle between the radius vector r*,^ and direction of the external 
magnetic field. 

In low magnetic field, when the Zeeman energy splitting uq = 'jHq, and DDI energy 
Ud = . / — J^^'i'i M are commensurable quantities, the time for reaching equilibrium state 

describing by Eq.(l) is determined by ^ which is of the order of r2 26 1. 



In high magnetic field, » Ud, the time t for reaching equilibrium becomes very long, 
t » T2 , due to the fact that the energy spectra of the Zeeman energy and the DDI energy 



are very different 



26|. 



We will consider the spin system placed in high magnetic field as two thermodynamic 
subsystems, the Zeeman and dipole-dipole ones, each of them is at thermodynamic equi- 



librium with own spin temperature and between which there is an energy transfer. This 
approach is based on the following reasoning. 

The total Hamiltonian can be presented as the sum of two terms each commuting with 
other , Hz and {[Hz, H^] = 0), and the term Hnd = H^d — H^, which does not commute 
with both Hz and H^. H^ is the secular part of dipolar Hamiltonian (3) 

Hd = -l^h J2 "^JkP^ (cos e^k) {Slzjlzk - Ijik) , (4) 

jk 

where P2 (cosOjk) = | (1 ~ 3cos^6'jfc), laj is the projection of the angular momentum op- 
erator of the j-th spin (j = 1, 2 , . . ., N) on the a - axis (a = x,y,z). At uq >> 
the effect of the dipolar interaction is small in comparison with the Zeeman interaction and, 
in the first-order perturbation, the secular part Hd (4) plays the main role in establishing 
equilibrium in each subsystem. This process is characterized by flip-flop terms in the secular 
part Hd , that ensures the density matrix in the following approximate form 

where (3z = and (3d = inverse Zeeman and dipolar temperatures, re- 

spectively, Z {f3z, (3d) = Tr {exp {—(3zHz — (3dHd)} is the partition function that encodes the 
statistical properties of the system in the quasi-equilibrium state. Separation of the spin 
system into the Zeeman and dipolar subsystems is consequence of the fact that the Zeeman 

/u\ /T\ d\nZ{(3z,^d) /„x 
{Hz) = uo (Iz) = ^ (6) 

and dipolar energies 

din Z{^,Pd) 

^""'^ = — — 

are separately constants of the motion , while both (Hz) and (Hd) depend on both thermo- 



dynamics parameters, (3z and (3d [26|. Results (5-7) have been conformed by experiments 

The time for reaching the quasi-equilibrium state described by Eq. (5) is of the order 
of T2. It is assumed that further relaxation process takes more time t > T2. This process 
consists in equalization of the temperature (3~^ and (3^^ and results in establishing of the 
equilibrium state which is characterized by density matrix (1) with Pz =(3d =(3- 

The approximate form of the density matrix does not include the non-secular term of the 
dipolar Hamiltonian, H^d ■ To take into account this term and to describe the relaxation at 



time t > T2, we will use below the method of the non-equilibrium statistical operator 



32| . 



III. ESTABLISHMENT OF EQUILIBRIUM 



During the time T2 the system loses the quantum coherence features and achieves a quasi- 
equihbrium state due to the flip-flop terms including in the Hamiltonian H^. For the spin 
system at low temperature, the decoherence time T2 typically ranges between nanoseconds 
seconds |33|. A typical valuem of T2 in a nuclear dipolar coupled spin system is 100 



ana 



34j | which is much shorter than the spin-lattice relaxation times of the Zeeman, T.^ and 



dipolar, Tid energies. Ti^ and Tid take values in the range from minutes to hours [26|]. We 
will study spin dynamics at T2 < t ^ min {Ti^, Ti^} when each subsystem can be considered 
as being at thermal equilibrium, and the spin-lattice relaxation process can be ignored. 
To derive an equation describing the establishment of equilibrium state (1), we take the 



density matrix in the general form 



32| 



Pne (t) = Zll exp ^ (^) - (^) f dre'^K, (r)^ j 



n = z,dd (8) 



where 

Zne (t) = Tr jexp ^ (t) if„ - /3„ (t) j' dre'^K,, (r)^ j | , (9) 
Kn (r) is the operator of the energy flux 

K„(r) = ^ (10) 

dr 

and transition to the limit e — > should be made after the calculation of the integral. 
According to the method of the non-equilibrium statistical operator , the operators under the 
integral in Eq. (10) are taken in the Heisenberg representation to be Hn (r) = e*^'^i/„e~*^^ 



32| 



Accounting that Yin (t) = 0, we rewrite density matrix (8) in the following form 32 1 



Pne = Z-^^l exp ^ /3fc (t) - (t) - /3, (t)) dre'^K (r) j . (11) 

Evolution of the spin system can be sufficiently completely described by the density 
matrix (11) with temporally depending ji^ and fid- To determine the temporal dependence 
of the inverse temperatures fi^ (t) and (3d (t) , let us calculate the average energy flux. This 
flux is given by averaging expression (10) using non-equilibrium density matrix (11) 
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(12) 

where (...) = Tr{pne---)- Eqs. (11) and (12) are valid at any temperatures of the sub- 
systems and can be apphed to consider dynamics of the dipolar coupled spin system and 
generation of the entangled state in the system. The quantities in Eq. (12), averaged with 
density matrix (11) , depend nolinearly on the inverse temperatures (3z and Pa and their 
difference /3z — /3d-Eq. (12) is solved using different approximations: the high-temperature 
approximation, <ti 1; the approximation of small values of the temperature difference, 
l/^z — <^ min or/and the approximation of small energy flux, \\K\\ <^ \\Hk\\ 

(||...|| denotes a norm of an operator) 32 1. 

As it is shown below, the entangled states can be achieved in the cases when a difference 
{f3d — 15 z) cannot be regarded as a small parameter. We cannot also use the high-temperature 
approximation and the expansion of (11) in powers of K . Hence, the perturbation theory 
methods cannot be applied effectively to study entanglement in the considered systems. 

We will solve numerically Eq. (12) for four-spin systems with various geometry: chain, 
circle, and rectangle. Examples of such systems are quasi-one-dimensional fluorine chains in 
calcium fiuorapatite Ca^F (PO^)^ [soj] , xenon tetrafiuoride with chemical formula XeF4^{its 
crystalline square planar structure was determined by both NMR spectroscopy and by neu- 
tron diffraction studies |40|, |4l|), and 1,4-dichlorobenzene CqH^CI with a rectangle planar 



structure. 

We use dimensionless units where the dipolar coupling constant of the nearest spins is 
chosen to be D = 1. The time unit is determined by the DDI energy between these spins 
and equals 100 fis . The ratio of Zeeman energy splitting to the dipolar energy is chosen to 
be — = 43 for all considered cases. We assume also that the angles Oj k = ^ for all pairs of 
spins in all cases. In our numerical calculation we used: a) for a chain all angles ipj^k = f 
and the coupling constants of spins the j-th and k-th are Dj^k = {j — k)^^; b) for a circle 
^12 = f , (P23 = = TT, V934 = ^, (fis = ^, ip24 = ^ and Dj^k = ^^^J^^ ! and c) for 

a rectangle = f , v^23 = '^u = tt, = ^, v^ia = arccos (^-^^, ^2a = arccos + tt, 
and D12 = -D34 = 1, D23 = = D13 = D24 = |. Fig. 1 shows evolution of the 
inverse temperatures (3z and /3d toward the equilibrium value of the four-spin systems in the 
form of a chain, a circle, and a rectangle. 
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These systems demonstrate qualitatively similar time dependences of j3z and j3d- At a 
high ratio ^ = 45 it is possible to mark out four stages in the j3d time dependences. These 
stages are clearly seen in the time dependence for a circle spin system (Fig. la, black 
line). The first stage, < t < 0.5 is characterized by a fast drop of the inverse temperature 
about twice . During the second stage, 0.5 < t < 2 , the rate of the [3d drop is decreased, 
and it is increased again at the next stage 2 < t < 2.5 . The dependence is substantivally 
non-exponential. Only at the last stage when of the temperature difference is small, the 
time dependence can be fitted by an exponential laws with the characteristic time close to 
r ~ 1. Of course, these stages can be marked out in the Pz time dependences. However 
the variations of the Zeeman temperatures are very small because the heat capacity of the 
Zeeman subsystem is high. The inverse temperature of the Zeeman subsystem j3z increases 
less than by 2.5% (Fig. lb), while the inverse temperature of the dipolar subsystem decreases 
70 times (Fig. la). 

With decreasing a ratio of the Zeeman energy splitting to the dipolar energy, the behave 
of the time dependence is changed. Only two stages can be marked out as a ratio ^ = 20 
(Fig. la, blue dash-dotted line). Both stages can be fitted by the exponential laws with 
characteristic times, r: r = 0.6 for the first stage < t < 0.5 and r = 3.3 the second one by 
. So, the decrease of the ratio of the Zeeman energy splitting to the dipolar energy leads 
to the increase of the characteristic time of achieving the equilibrium state. 



IV. ENTANGLEMENT AT A LOW ZEEMAN AND DIPOLAR TEMPERATURES 

Several parameters have been proposed for characterizing the entangled state of a spin 
system: the von Neumann entropy, entanglement of formation, log negativity, concurrence 



3ntrop\ 



of a pair of spins, and etc. [2|, |3|, |35H38| . We will characterize the entangled states by the 
concurrence between two, m-th and n-th, spins which is defined as jssi 

Cmn iPz, Pd) = max {qmn {f^z, Pd) , 0} , (14) 

with {/3z, Pd) = >!^nl (/?., Pd) - (/?., (5d) - \^^n (/?., Pd) - X[l (/?„ (3d). Here A^^l (/?„ /?,) 
{k = 1,2,3,4) are the square roots of eigenvalues, in the descending order, of the following 
non-Hermitian matrix: 

■^mn i,(3zj Pd) Pmn (^(^zj f^d 



where pmn {f^z, Pd) is the reduced density matrix. For the m-th and n-th spins, the reduced 
density matrix pmn is defined as pmn = Trmn (p) where Trmn (•••) denotes the trace over 
the degrees of freedom for all spins except the m-th and n-th spins. In Eq. (15) pmn is 
the complex conjugation of the reduced density matrix p^n and cXy is the Pauli matrix 
/O -^\ 

Cly = \ • For maximally entangled states, the concurrence is C^n = 1 while for 

V V 

separable states Cmn = 0. 

Figs. 2 and 3 present the results of numerical calculations of the concurrence Cmn in the 
four-spin circle (Figs. 2a and 3a), rectangle (Figs. 2b and 3b), and chain (Figs. 2c and 3c) 
. Dependences of concurrencies on inverse temperatures (3z and /^d at ^ = 45 are presented 
in Figs. 2. Here the dashed red lines present concurrencies in the equilibrium states ( 
f3z= l^d) and the dotted yellow lines show dynamic concurrency variations at achieving the 
equilibrium state. 

A. Conclusion 

We have investigated the entanglement evolution during establishing an equilibrium in 
four-spin dipolar coupled systems with various structure: circle, rectangle, and chain. The 
concurrence was considered as a measure of entanglement. Entanglement is observed in non- 
equilibrium states ( 7^ /?d) and can appear when both temperatures much higher than 
the temperature of the entanglement appearance in the equilibrium state. This confirms 
the results of previous works that non-equilibrium conditions are more favorable for the 
generation of entangled states as equilibrium states 



23 



24| . Entangled states appear and 



are stored at lower dipolar temperatures relative to the Zeeman one( < /S^^). We found 
that the concurrence non-monotonically depends on both temperatures (Fig. 2). 

Considering the establishment of the equilibrium state, we showed that the well-known 
high-temperature approximation is not applicable for a homonuclear spin system at low tem- 
perature, and the temporal dependence of the inverse temperatures of the spin subsystems 
is given by a solution of non-linear equations. 

We demonstrated that the concurrence dependences on time is non-monotonic during the 
equilibrium state establishment process (Fig. 3) and entanglement fades long before equilib- 
rium is established (compare Figs. 1 and 3). The behavior of entanglement depends strongly 
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on the ratio of the Zeeman energy to the dipolar energy. At high ratio, the concurrence in 
the system decreases significantly till the residual very small quantity in time of the order 
of T2. At low ratio, the establishment of equilibrium (Fig. la and b, blue dash-dotted fine) 

and the decay of entanglement (Fig. 3a, green dotted line) take time up to 1 ms. 

Obtained results open an effective way to the experimental testing of entanglement in 
spin systems by using non-equilibrium states to generate the entangled state at higher tem- 
peratures. 
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Figure Captions 

Fig. 1 Equalization of the inverse spin temperatures of the dipolar (a) and Zeeman (b) 
subsystems in a four spin circle (solid black line), chain (dashed red line), and rectangle 
(dotted yellow line) at ^ = 45. Blue dash-dotted hne - the inverse spin temperatures of 
the dipolar (a) and Zeeman (b) subsystems in a four spin circle at ^ = 20 . 

Fig. 2 The concurrence in four spin system as a function of the inverse temperatures /3z 
and /Sd- (a) Cu in the circle, (b) C13 in the rectangle , and (c) C14 in the chain. Dashed 
red line - equilibrium state with I3z= Pd- Dotted yellow line - evolution of the concurrence 
under equalization of the Zeeman 13^ and j3d inverse temperatures due the influence of the 
non-secular term Hnd- 

Fig. 3 Time dependence of concurrence Cmn in ^ four spin system at ^ = 45: (a) circle, 
(b) rectangle, (c) chain. Sohd black hne - C12, dashed red hne - C13, and dotted blue hne 
- C14. Green dotted line - C12 in a four circle at ^ = 20. 
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